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Abstract 

Motivated by recent progress of theory and experiment on the anomalous magnetic mo- 
ment of the electron, a e , we update the hadronic contributions to a e . Using our up-to-date 
compilation of e + e~ — > hadrons data, we find the leading order hadronic contribution 
a had, lo, VP _ ^ ggg _|_ o.010 exp ± 0.005 ra d) • 1CT 12 and the next-to-leading order hadronic 
contribution a ^d,NLO,VP = (_ .2234 ± 0.0012 exp ± 0.0007 rad ) • lCT 12 , where the first and 
second errors are from the error of the experimental data and the uncertainty in the treat- 
ment of radiative corrections, respectively. These values are compatible with earlier evalu- 
ations by other groups, but have significantly improved uncertainties due to the more pre- 
cise input data used. We also update the leading order hadronic contribution to the ground 
state hyperfine splitting of muonium, obtaining Ai/^u = (232.68±1.25 exp ±0.72 rad )Hz. 
This value is consistent with the most precise evaluation in the literature and reduces its 
error by a factor of two. 

1 Introduction 

With the start of the LHC, the hunt for physics beyond the Standard Model (SM) at the 
energy frontier has entered a new era, though so far no direct signal for 'new physics' has been 
observed. At the same time, experiments at lower energies are measuring SM parameters with 
unprecedented accuracy and are becoming ever more sensitive to quantum effects, possibly from 



physics beyond the SM. A prime example for this is the anomalous magnetic moment of leptons 
PQ, caused by loop effects mainly from Quantum Electro-Dynamics (QED) but also influenced 
by the strong and weak force and possibly by effects from new physics. For the muon, the 
anomaly a M is sensitive to all sectors of the SM and has been measured with 0.5 ppm accuracy, 
using a muon storage ring at Brookhaven [2]. This measurement, when compared to the SM 
prediction of a^, shows a discrepancy of 3.3 a [3J HI IS] which may be a sign for new physics^ 

The electron's anomalous magnetic moment, a e , has been measured at Harvard via atomic 
spectroscopy [TJ, 



This corresponds to a 0.24 ppb accuracy, hence a e is known even much more precisely than 
a M . However, due to the small mass of the electron, a e is less sensitive to quantum effects 
from higher mass scales and is mainly generated by QED effects, though it also receives small 
contributions from the strong and weak interaction. It is therefore an ideal place to determine 
the fine structure constant a. 

Another precision observable is the ground state hyperfine splitting (HFS) of muonium. 
Similar to a e , it is sensitive to quantum effects at low mass scales and mainly QED dominated. 
It is very useful to determine the electron-to-muon mass ratio and hence the muon mass. 

To determine SM parameters or to become sensitive to new physics, it is crucial that the 
theoretical predictions of such precision observables have accuracies similar to or smaller than 
the experimental uncertainties. For the lepton anomalies and also the HFS of muonium this 
has been achieved by higher-order calculations in all sectors of the SM (and various extensions 
of it) [31 [8j [9]. For a e , the group of Kinoshita and collaborators have, after a many- year effort, 
completed calculations including full five-loop effects [TQJ [TTJ [131 EEU EH EES] • These results 
have reduced the error of the theoretical prediction of a e by more than four-fold. It is therefore 
timely to provide an up-to-date evaluation of the hadronic contributions. In this short letter 
we calculate the hadronic contributions due to vacuum polarisation effects for the anomalous 
magnetic moment of the electron and also for the hyperfine splitting of the ground state of 
muonium. 

2 Anomalous magnetic moment of the electron 

At present, the most precise determination of the fine structure constant is through the measure- 
ment of the anomalous magnetic moment of the electron, a e , using the Harvard measurement 
[7J quoted in Eq. flU) above. With their new complete five-loop and updated and improved four- 
loop QED contributions, and using results from [TTJ [THJ [T^J [20] for the hadronic and electroweak 
(EW) contributions, Kinoshita et al. arrive at [TO] 



lr rhe tau lepton's anomaly, a T , is even more sensitive to physics at higher energy scales than a^, but is very 
difficult to measure with high precision due to the short lifetime of the tau. See [5] for a recent review and SM 
prediction of a T . 
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corresponding to a 0.25 ppb accuracy. This value is slightly different but compatible (at the 
level of 1.3 a) with the second best determination of a which is using Rubidium atoms [2H 122] 
and reads cr 1 (R.b) = 137.035 999 049 (90) [0.66 ppb]| The four errors quoted in ([2]) stem from 
the four- and five-loop QED contributions, from the hadronic and electroweak effects, and from 
the uncertainty of the experimental value for a e . The group of Gabrielse is working on fur- 
ther improvements of the measurement of a e and also a similar determination of the positron's 
anomaly which, assuming CPT invariance, can be combined with the electron's anomaly similar 
to the case of the muon. However, currently the experimental error of a e is dominating the 
theoretical uncertainties, which have been reduced dramatically by the new QED calculations 
of Kinoshita and collaborators. Their four- and five-loop errors, i.e. the first and second error in 
(J2J), are mainly coming from the statistical uncertainties of the numerical Monte Carlo integra- 
tions and can be improved with increased computational efforts when needed. To understand 
the composition of the third error, let us list the hadronic and electroweak contributions as 
used in [10J: 

a had,LO,VP = i. 8 75 (is). iQ- 12 , (3) 
a had,NLO,VP = _ .225 (5) ■ KT 12 , (4) 

a had,l-by-l = Q 035 ( 1Q ) . 10 -1 2) ( 5 ) 

af w = 0.0297 (5) -10- 12 . (6) 

Here, similar to the case of the anomalous magnetic moment of the muon, the error from the 
electroweak effects is negligible compared to the uncertainties from the hadronic contributions. 
The latter are dominated by the leading order (LO) and next-to-leading order (NLO) hadronic 
vacuum polarisation (VP) corrections. For a e , the contributions from light-by-light (1-by-l) 
scattering diagrams are rather small, but they still add significantly to the uncertainty of the 
Og ad . A further scrutiny of the model calculations, and possibly also 'first-principle' determina- 
tions based on lattice calculations, are under way for the anomalous magnetic moment of the 
muon. This in turn will also allow to improve the estimates for a k ad > 1_b y~i i n the following, 
we will discuss our up-to-date determination of a h ad > LO > VP an d a h ad > NLO > VP . Note that while 
the hadronic cross section data used as input in these calculations have improved significantly 
in recent years due to both direct scan measurements and via the method of radiative return 
(see e.g. [21] for a recent review of the field), no up-to-date determination of the VP induced 
contributions is available so far. 

2.1 Leading order hadronic VP contributions 

The calculation of the leading order hadronic VP contributions to a e is very similar to the case 
of the muon and uses a dispersion integral over the hadronic cross section times a well-known 
kernel function K^(s), see e.g. the detailed discussion in [25J: 

i r°° 

4 ad ' L °' VP = ^3 / 2 & K t (s)o°U*) , (7) 



2 Note that this value is a big improvement from the previous measurement by the same group 
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Energy region, remarks 


contribution to a e 


contribution to Ai/j^' P 


Zm^ — U.oUO Lrev ^Lllr 1 , ZTT) 

3m n - 0.66 GeV (ChPT, 3vr) 
- 0.60 GeV (ChPT, vr 7 ) 
m v - 0.69 GeV (ChPT, 777) 
— > unaccounted modes 


U.UUol ± U.UUU1 


0.0004 ±0.0000 


0.0001 ±0.0001 


U.zo it U.U1 


0.04 ±0.00 


0.01 ±0.01 


0.305 — 1.43 GeV (exclusive channels) 

1.43 — 2 GeV (excl. data + isospin relations) 

2 — 2.6 GeV (inclusive data) 

z.o — 0.10 kjiQ v ^p^LyU wren rjri/oii error j 

S 7S — 1 1 DQ OpV finrl data"! 


1.0000 ± U.UUyz 

0.0846 ± 0.0027 
n 0^78 -i- n nm k 

u.uzo( it u.uuuy 
n oqac. 4- n nnn^ 

\J .KJO'-kO _I_ U.UUL/tJ 


190.10 ± l.UY 
13.41 ±0.43 

0.04 it U.ZO 

a qa _i_ n 1 7 

4.o4 ± U.li 

7 47 4. n 1 n 

I .4: I _I_ U.1U 


J/V^ + ^ 
T(1S-6S) 


0.0185 ±0.0004 
0.0002 ± 0.0000 


3.53 ±0.07 
0.06 ±0.00 


11.09 -00 GeV (pQCD) 


0.0049 ± 0.0000 


1.37 ±0.00 


Total 


1.866 ±0.010exp 


232.68 ± 1.25exp 



Table 1: Contributions to a had,LO,VP un ^ s Q f 10~ 12 ) and to Az/ M a u ' (in Hz) from different 
energy regions, obtained with the data compilation as used in [3]. The first four lines give our 
predictions of contributions close to threshold where no data are available and are based on 
chiral perturbation theory (ChPT), see [25] for details. For 2.6 < yfs < 3.73 GeV perturbative 
QCD (pQCD) with errors comparable to those of the latest BES data [26] in this energy region 
is used. In the region below 2 GeV the sum of exclusive channels, supplemented by isospin 
relations for channels where no data are available, is used, see [3] for details. 



where I = e, /j, and cr° ad (s) is the undressed (i.e. excluding VP corrections) total hadronic cross 
section. If we define the function Ke(s) by 

k t {a) = — 2 K t {s) , (8) 

m i 

then Ki(s) is a monotonically increasing function of order one with Ki(s) — > 1 for s — > 00. In 
the case of a e , the kernel K e (s) is very close to one throughout. In fact, K e (s) can be expanded 
in terms of m? e /s as 

4(^i + ( 3 ,4 + H5)^ + (^.„^), (9 ) 

and the deviation of K e (s) from one is almost always negligible. (In the case of a M , K^{s = 
77^.0) = 0.40, Kfj,(s = 4:m^ ± ) = 0.63, and K^s) — > 1 for s — » 00.) Hence the low energy 
contributions to the dispersion integral are even more important than in the case of the muon. 
Using our up-to-date comprehensive compilation of hadronic cross section data [1] we obtain 
contributions to a e from different energy regions as displayed in Tab. [TJ The errors given in 
Tab. [T] contain the statistical and systematic uncertainties from the input data including corre- 
lations over different energies for the various hadronic final states which are added incoherently. 
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In addition, we have to take into account the uncertainties in the treatment of radiative cor- 
rections from final state radiation and VP effects, see [U [25] for a detailed discussion. These 
are conservatively estimated to be 

5al s R = 0.005 , 5aY = 0.002 (10) 

and have to be added in quadrature to the error of the total Og ad ' LO ' VP which hence reads 

a had,LO,VP = ( Lg66 ± Q11 ) . 1Q -12 _ ^ 

This number is to be compared with the result from [17] given in Eq. ([3]), which has been 
used in [TO] and [22]- As expected, the mean value has changed only slightly, but the total 
error is significantly improved by about 40% due to the improved hadronic cross section data. 
Compared with the more recent evaluation given in [8], which is a ^ ad > LO > VP = (1.860 ± 0.015) • 
10~ 12 , the agreement is very good with a further but less significant improvement in the error. 



2.2 Next-to-leading order hadronic VP contributions and total 

Contrary to the next-to-leading order (NLO) contributions induced by so-called hadronic light- 
by-light scattering diagrams, the NLO corrections including hadronic VP diagrams can be 
calculated with help of dispersion integrals, see e.g. (H [HI [25]. Due to the small mass of the 
electron, diagrams with two hadronic VP insertions or one hadronic and one VP insertion from 
a heavy lepton (muon or tau) are strongly suppressed, and the diagrams with one hadronic VP 
insertion and one additional photon or electron loop are practically the only relevant contribu- 
tions. With our latest compilation of hadronic cross section data we obtain the value 

a had,NLO,VP = (_ .2234 ± 0.0012 cxp ± 0.0007 md ) • 10~ 12 , (12) 

where the uncertainty due to the statistical and systematic errors of the hadronic cross section 
data used as input and the additional error due to radiative corrections applied to the data 
are given separately. This result is close to the result quoted in Eq. (TjJ [IS] , but has a much 
smaller error. These NLO corrections lead to a reduction of the hadronic LO VP corrections 
by about 12%, so should not be neglected. 

Combining the results from Eqs. (flT!) and (fl2l N with a^ ad ' 1-137-1 from ([5]) \19^ , our estimate 
of the hadronic contributions to the electron's anomaly is 

a^ d = (1.678 ± 0.014) • 10~ 12 . (13) 

Equations ffTTl) and ffl2l are the main result of this letter. Given the current uncertainties of 
the QED result for a e and the uncertainty of the experimental measurement, the improvements 
in the hadronic contributions may not look dramatic and do not affect the determination of a 
significantly. However, anticipating future improvements for a^ ED and ag Xp , it is important that 
the hadronic contributions are now much better under control and will not limit the indirect 
determination of a for the foreseeable future. 

3 The LO and NLO VP contributions are, to a very good approximation, totally correlated, which is taken 
into account accordingly in the error combination. 

4 This value for the light-by-light contributions is compatible with the evaluation by Jegerlehner and Nyf- 
feler [SJ|27], who quote a^ 1 '^' 1 = 0.039 (13) • lO" 12 . 
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3 Hyperfine splitting of muonium 



Another high-precision observable is the HFS of the ground state of muonium. Comparing 
its measured and predicted values is currently the best method to determine the electron-to- 
muon mass ratio and hence the muon mass [22] ■ The HFS is mainly QED dominated, but 
also receives higher-order contributions due to hadronic and EW interactions, see 1221 [2H] for 
detailed reviews and further references. 

The experimental value of the muonium HFS is [29, 30J: 



z/ Mu (exp) = 4 463 302 776(51) Hz , (14) 

while the theoretical value quoted in the 2010 version of CODATA [22] is 

z/ Mu (theory) = 4 463 302 891(272) Hz , (15) 

where the uncertainty is dominated by that of the mass ratio m e /m^. According to Ref. [22J, 
among the theoretical uncertainty of 272 Hz, the uncertainty of the QED contributions is 98 
Hz and that of the hadronic contribution 4 Hz. The hadronic contribution quoted in Ref. [22J 
is 

Ai£? = 236(4) Hz, (16) 

which is the sum of the hadronic VP contribution Ai/ Mu ,VP = 231.2(2.9) Hz [31] and the 
hadronic higher-order contribution Az/ Mu ' H = 5(2) Hz [31] . Contrary to the case of the lep- 
ton anomalies, for the muonium HFS the hadronic light-by- light contributions are completely 
negligible at the current level of accuracy [32] . While the error in the theoretical prediction of 
the muonium HFS is dominated by the estimates of unknown higher-order QED contributions, 
the hadronic contribution is not negligible and the error quoted in ( 1161) is just about one order 
of magnitude below the current experimental error. In fact, there is a planned experiment to 
measure the muonium HFS at J- PARC [33], which aims at reducing the experimental uncer- 
tainty by a factor of two or more compared to Eq. (fT4T) FI We therefore take the opportunity 
to improve the hadronic contributions. The hadronic VP contributions to the muonium HFS, 
Ay Mu |VP , have previously been evaluated by a number of groups [HH ESI ESI EH EH1 ESI l4"0"] . 
They can be written SIS db dispersion integral 

A had. VP 1 m < 



POO 

-v F / ds K Mu (s)0-£ ad (s), (17) 



where vp is the so-called Fermi energy, 



16 2 m, 

vf = —Roc® — 
3 m 



1 + 



m 



l -3 



nip 



(18) 



5 Also note that in their most recent work on QED corrections to the HFS of muonium, Eides and Shelyuto 
quote an uncertainty of about 10 Hz as current goal for the theoretical uncertainty of the HFS |34) . 
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where is the Rydberg constant, = 3 289 841 960 364(17) kHz [22]. The explicit form of 
the kernel function K Mu (s) is given e.g. in Ref. [39]. After correcting typos in Ref. [39] it reads 



AW*) H " 2 + 2 ) j + + f ) hl * - *■ (fOT ' *= 4, " y ' (19) 

" (*% + 2 ) /31n (£§) + (4% + l) - 1, (for s > 4mJ) , 



where 



^s-^'t^-i^Eg, (20) 

r r V s 

The kernel ifu u (s) is a monotonically decreasing function with Km u (s = w 2 ) = 5.53 and 
-K"mu(<s = 4m 2 ±) = 1.98, approaching zero as s — » 00. For large s, an expansion in terms of 



m^/s is useful: 



9 , ml 15 \ m?. „„ / mf. , mf 



Klb{t) = [ U-f + T \-f + Ol^-f] . (21) 



Basically the dependence of Kuu{s) on s is similar to that of K^(s). Since the leading term in 
the above expansion is irrvi/s) ln(s/m 2 ) rather than (m^/s), it puts only slightly more emphasis 
on higher energies compared to the hadronic VP corrections to a M . 

Using the same input as above, we obtain the hadronic contributions from different energy 
regions as listed in Tab. [TJ As in the case of a e above, the error displayed in the last line of 
Tab. [T] contains only the uncertainties of the experimental data. Adding an error from the 
conservatively estimated uncertainty due to radiative corrections, ±0.72 ra d Hz, our final result 
for the hadronic VP contributions to the muonium HFS reads 

Ai^f- VP = (232.68 ± 1.44) Hz . (22) 

Our result is compatible with and, as expected, considerably more accurate than the previous 
result quoted in Eq. ffl6|) . Note however, that so far no attempt has been made at a complete 
calculation of higher-order hadronic VP corrections which have been estimated to be of the 
same order or bigger than the error of the leading order ones 



4 Conclusions 

We have used our comprehensive compilation of hadronic cross section data to determine the 
hadronic vacuum polarisation contributions to the anomalous magnetic moment of the electron, 
a e , to leading and next-to- leading order. We have also evaluated the hadronic VP contribu- 
tions to the hyperfine splitting of the ground state of muonium. Our main results are given in 
Eqs. (|lip . ( 1T2|) and f l2"2"j) . While the changes of the central values compared to earlier determi- 
nations are modest, the corresponding uncertainties have been significantly improved and will, 
for the foreseeable future, not affect precision determinations of physical constants from these 
observables. 
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